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1 Introduction 



One of the challenges in quantum chromodynamics (QCD) is the understanding of 
hadronic high-energy reactions. Asymptotic freedom allows perturbative calcula- 
tions of hadron-hadron interactions at large momentum transfer. Computer simula- 
tions on Euclidean lattices give access to non-perturbative QCD and allow to inves- 
tigate static properties of hadrons such as the hadron spectrum and confinement. 
Less is known about the non-perturbative structure of hadron-hadron scattering as 
non-perturbative methods working in Minkowski space-time are required. 

In this work, we investigate the perturbative and non-perturbative QCD struc- 
ture of the dipole-dipole scattering amplitude in momentum space. We use the 
recently presented loop-loop correlation model (LLCM) |I| which combines pertur- 
bative and non-perturbative QCD contributions and allows to compute cross sections 
of hadron-hadron, photon-hadron and photon-photon reactions in a unified descrip- 
tion. The model is based on a functional integral approach developed for parton- 
parton scattering in the eikonal approximation f^,|| and extended to gauge-invariant 
loop-loop scattering The S"-matrix factorizes into the universal correlation of 

two light-like Wegner- Wilson loops and reaction-specific wave functions. The light- 
like Wegner- Wilson loops describe color-dipoles given by the quark and antiquark 
in the meson or photon and in a simplified picture by a quark and diquark in the 
baryon. The size and orientation of the color-dipoles in the hadrons and photons are 
determined by appropriate light-cone wave functions. The perturbative correlation 
is described by lowest order gluon exchange and the non-perturbative correlation by 
the stochastic vacuum model (SVM) 0. The unitarization of both contributions 
obtained in the approach of Berger and Nachtmann || is important for ultra-high 
energies and non-forward scattering |T[ and is not discussed here. 

The SVM gives confinement for non-Abelian gauge theories due to flux-tube 
formation of color fields between the quark and antiquark of a dipole PHTT|. This 
flux-tube or string enters via the essential SVM assumption that non-perturbative 
interactions are most adequately described in terms of correlators of gluon field 
strengths instead of gluon potentials as in perturbative QCD. Therefore, the line 
integrals over the gluon potentials in the Wegner- Wilson loops of the functional in- 
tegral approach are transformed into surface integrals over gluon field strengths 
using the non-Abelian Stokes theorem |3],[T2|]. In line with our recent work |]J we 
use minimal surfaces, planar surfaces bounded by Wegner- Wilson loops, which are 
usually used to obtain Wilson's area law in Euclidean space-time |f7|,|II|. Mini- 
mal surfaces allow us to show for the first time the non-perturbative structure of 
the dipole-dipole scattering amplitude in momentum space. The previous pyramid 
mantle choice for the surfaces p|, |3|, |b ^lTfl did not allow such an analysis. 



2 



The perturbative contribution to dipole-dipole scattering in our approach is the 
known two-gluon exchange [H|[2(| between the two dipoles. The non-perturbative 



contribution comes out as a sum of two parts: The first part describes the non- 
perturbative interaction between the quarks and antiquarks of the two dipoles and 
exhibits the same dipole factors as the perturbative contribution. The second part 
represents the interaction between the strings of the two dipoles. It shows a new 
structure different from the perturbative two-gluon exchange. Nevertheless, we 
would like to empasize that \k± \ - factorization is valid also for the non-perturbative 
contribution. 

The string confining the quark- antiquark dipole of length \f D \ has an integral 
representation which sums stringless dipoles of sizes £\r D \ with < £ < 1 and dipole 
number density n(£) = l/£ 2 - Consequently, the string-hadron scattering reduces 
to an incoherent superposition of dipole-hadron scattering processes. This decom- 
position of the string into stringless dipoles allows us to extract the microscopic 
structure of the unintegrated gluon distribution of hadrons and photons, J-h(x, kj_), 
from our dipole-hadron and dipole-photon cross section via \k±\ -factorization. 

The unintegrated gluon distribution of hadrons and photons J r f l ( basic, 
universal quantity convenient for the computation of many scattering observables 
at small x. It is crucial to describe processes in which transverse momenta are 
explicitly exposed such as dijet |21| or vector meson [pj| production at HERA. 



Its explicit | k± \ - dependence is particularly suited to study the interplay between 
soft and hard physics. Moreover, the unintegrated gluon distribution is the central 
object in the BFKL [E3] and CCFM evolution equations. Upon integration over 



the transverse gluon momentum \k±\ it leads to the conventional gluon distribution 
xGh{x,Q 2 ) used in the DGLAP evolution equation [EH]. 



The outline of the paper is as follows: In Sec. we review the loop-loop correla- 
tion model U and give the model parameters. In Sec. [| we investigate the pertur- 
bative and non-perturbative structure of the dipole-dipole scattering amplitude in 
momentum space. In Sec. f|, we present the decomposition of the string into dipoles, 
elaborate the total dipole-hadron cross section, and extract the unintegrated gluon 
distribution. Numerical results for the unintegrated gluon distribution of protons, 
pions, kaons, and photons are given in Sec. [5] together with the integrated gluon 
distribution of the proton xG p (x, Q 2 ). In Sec. we compare our results for the un- 
integrated gluon distribution of the proton with those obtained in other approaches. 
We close with a summary of our results. 
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2 The Loop-Loop Correlation Model 



Recently, we have presented a loop-loop correlation model to compute high-energy 
hadron-hadron, photon-hadron, and photon-photon reactions JIJ . Based on the func- 
tional integral approach to high-energy scattering in the eikonal approximation |2|-|6| . 
the T-matrix element for the elastic scattering of two particles 1 and 2 at transverse 
momentum transfer cjj_ (t = —q 2 ) and cm. energy squared s reads 

T(s,t) = 2is Jd 2 b ± e^ Jdz l d 2 r l Jdz 2 d 2 r 2 

x |^i Oi,^)! 2 \ip 2 {z 2 , r 2 )\ 2 1 - S DD (b ± ,z 1 ,r 1 ,z 2 ,f 2 ) , (2.1) 

where the correlation of two light-like Wegner- Wilson loops, the loop-loop correlation 
function^ 

Sdd^zi,?!,^) = (W[C 1 ]W[C 2 ^ G , (2.2) 

describes the elastic scattering of two color-dipoles (DD). The path of each color- 
dipole is represented by a light-like Wegner- Wilson loop PB| in the fundamental 
representation of SU(3) 



^[Cy^TrPexp 



ig <p dz^Q^z) 



(2.3) 



where Tr is the trace in color space, g the strong coupling, and G^(z) = Q®(z)t a 
the gluon field with the SU (3) group generators t a that demand the path ordering 
indicated by V . Physically, the Wegner- Wilson loops (|2.3|) represent the phase 
that quarks and antiquarks acquire along the light-like trajectories Cj in the gluon 
background field. The color-dipoles have transverse size and orientation r*j. The 
longitudinal momentum fraction of dipole i carried by the quark is Zj. The impact 
parameter between the dipoles is [[ljj 



b± = f lq + (1 - zi)fi - r 2q - (1 - z 2 )f 2 = r lcm - r 2cm , (2.4) 
where r iq (rig) is the transverse position of the quark (antiquark), fi = r iq - — r iq , and 



1 i cm « 6 i ' iq 



+ (1 — Zi)fi q is the center of light-cone momenta. Figure [TJ illustrates 
the (a) space-time and (b) transverse arrangement of the loops. The QCD vacuum 
expectation value (. . .)g in the loop-loop correlation function ( |2.2j ) represents func- 
tional integrals in which the functional integration over the fermion fields has 
already been carried out as indicated by the subscript G. The model we use for the 



5 More generally, S DD = {W[Ci]W[C 2 ])gI ({W[C{\) g {W[C 2 })g)- For light-like Wegner- Wilson 
loops, however, one obtains {W[Cia])a = 1 as shown explicitly in 0. 
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Figure 1: High-energy dipole-dipole scattering in the eikonal approximation represented 
by Wegner- Wilson loops: (a) space-time and (b) transverse arrangement of the Wegner- 
Wilson loops. The shaded areas represent the strings extending from the quark to the 
antiquark path in each color dipole. The thin tube allows to compare the field strengths 
in surface S± with the field strengths in surface S^- The impact parameter b± connects 
the centers of light-cone momenta of the dipoles. 
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QCD vacuum describes only gluon dynamics and, thus, implies the quenched ap- 
proximation that does not allow string breaking through dynamical quark- ant iquark 
production. 



The fi and Z{ distribution of the color-dipoles is given by the wave functions ipi 
that characterize the interacting particles. In this framework, the color-dipoles are 
given by the quark and antiquark in the meson or photon and in a simplified picture 
by a quark and diquark in the baryon. We use for hadrons the phenomenological 
Gaussian wave function fl8|,[27| and for photons the perturbatively derived wave 
functions with running quark masses rrif(Q 2 ) to account for the non-perturbative 
region of low photon virtuality Q 2 [15] as discussed explicitly in Appendix |A]. 



The computation of the loop-loop correlation function fl2.2|) in the matrix cu- 
mulant approach of Berger and Nachtmann [R) has been reviewed in detail in 
The main steeps of this computation are the transformation of the line integrals into 
surface integrals with the non-Abelian Stokes theorem, a matrix cumulant expan- 
sion, and the Gaussian approximation of the functional integrals. In line with other 
two-component (soft + hard) models |T6|. [T7|, |28| pOfl an d the different hadronization 
mechanisms in soft and hard collisions, we have further demanded that the pertur- 
bative and non-perturbative contributions do not mix on the amplitude level 
The resulting T-matrix element at the reference cm. energy ^/sq reads 



T(s ,t) = 2is d 2 b ± e^ d Zl d 2 n dz 2 d 2 r 2 \M^n)\ 2 \M z 2,r 2 )\ 



x 



,2 1 P \ I NP \ 1 2 p \ 2 NP 

1 COS — Y COS — Y COS — Y COS — Y 

3 \3 A / V3 A / 3 V3 A / V3 A 



(2.5) 



where the arguments of the functions x F,NF iP^i z ii n ; z 2, ^2) are suppressed to lighten 
notation. The x~ functions defined as 



X 



,P,NP 



7T 



[ da^( Xl ) [ da^(x 2 )F^f(x 1 ,x 2 ,o;C xlo ,C X20 ) (2.6) 
JSi Js 2 



represent integrals over the minimal surfaces Si <2 that are built from the areas 
spanned by the quark and antiquark paths Ci >2 of the interacting dipoles and the 
infinitesimally thin tube which connects the surfaces S\ and S2 as shown in Fig. |l|. 
The integrand Fj^^f is the gauge-invariant bilocal gluon field strength correlator 



4vr 2 



0JL(°j X U Cx l0 )Gpa(0, x 2 ; C X20 ) 



P,NP 



G 



Ffivpa ( x i5 x 2, o; C X1Q , C X20 ) 



(2.7) 
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where the gluon field strengths Qp, v {xi) are parallel transported to the common 
reference point o along the path C Xi0 



G^(o, Xi] C Xt0 ) = o; C Xi0 ) 1 G fMl ,(xi)®(x, o; C XiQ ) 
with the QCD Schwinger string 



(2i 



o; C x . ) = Pexp 



CxaO 



(2.9) 



to ensure gauge invariance in the model. In (|2~7l ), gives the perturbative 

(P) physics (short-range correlations) described by perturbative gluon exchange and 
Fuvpa ^ ne non-perturbative (NP) physics (long-range correlations) modelled by the 
stochastic vacuum model (SVM) J7|. 

In leading order in the strong coupling g, the perturbative bilocal gluon field 
strength correlator F pvpa is gauge-invariant already without the parallel transport 
to a common reference point and depends only on the difference z := X\ — X2, 



F p 



fiupcr 



7l 2 2 



d , _d_ 

z u9fip Zp9p,cr) T \ Z p9v 



dz, 



z«g, 



up J 



Dp(z 2 



(2.10) 



^ 2 r d ^k 



-ikz 



tt 2 J (2vr) 4 



kuk a 9p,p k v k p g pu + kpk p g va kpk a g up D p (k 



Here, we have introduced an effective gluon mass mg = m p = 0.77 GeV as infrared 
regulator in the perturbative correlation function 



d 



D' P (k 2 ) :=— I d*zD P (z 2 )e 



„2\ „ikz 



k 2 — 171q 



(2.11) 



and a parameter M 2 = 1.04 GeV 2 which freezes the running coupling^ in the 
quenched approximation at the value a s {k"j_ = 0) = 0.4 [11], 



a s (k\ 



9\kl) 



4vr ll\n[(k 2 +M 2 )/A 2 QCD ] ■ 



(2.12) 



If the path connecting the points X\ and %2 is a straight line, the non-perturbative 
correlator F pi f pa depends also only on the difference z := x\— x 2 - Then, the most gen- 
eral form of the correlator in four-dimensional Minkowski space-time that respects 

6 Only the running coupling as a function of the transverse momentum is needed in the final 
result. 
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translational, Lorentz, and parity invariance reads @-§[ 



FZM ■= fZ ( « c \z) + FZ¥(z) (2-13) 

F^ c) {z) = 6( ^ 2 Z~Y) idz' ( z °9w~ z p9^ + q^- ( z p9ua-z a g vp )) D^z 2 ) (2.14) 

G2 (1 — K) f dk —ikzfi 1 11 11 ii \ n//72\ 

= ~~ 6(iV 2 lj J (27r) 4 6 y^f^^d^p ~ kvkp9iia + kpKpQua ~~ kpkoQvp J ^Al™ ) 

F^pV{z) = 2 _ [g^gua - g^gup) D{z ) 

^ " Gtofc* - V^p) / 7^1 ^(^) (2-15) 



3(iV c 2 -l) J (2tt) 

and was originally constructed in Euclidean space-time [0]. In all previous appli- 
cations of the SVM, this form depending only on xi and has been used. New 
lattice results on the path dependence of the correlator [31| show a dominance of the 
shortest path. This result is effectively incorporated in the model since the straight 
paths dominate in the average over all paths. 

In ( |2.13|) , D\ and D are the non-perturbative correlation functions in four- 
dimensional Minkowski space-time and 

D[(k 2 ) -=^j d^zD^z 2 ) e ikz . (2.16) 

The Minkowskian correlation functions D\ and D are derived |l| from the simple 
exponential correlation functions in four-dimensional Euclidean space-time 

Df (Z 2 ) = D E (Z 2 ) = exp(-\Z\/a) , (2.17) 

that are motivated by lattice QCD measurements of the gluon field strength corre- 
lator F^f pa (Z) |32],|33|]. The parameters are the correlation length a = 0.302 fm, the 

gluon condensate G 2 := (£sG^(°)^( )) = 0.074 GeV 4 , and the relative weight of 
the two different tensor structures k = 0.74. 

The correlator (|2.13|) is a sum of the two different tensor structures F^J}^ (z) 
and Fp, V pa(z). The first ( |2.14D is characteristic for Abelian gauge theories and 



does not lead to confinement, while the second ( |2.15| ) can only occur in non- Abelian 
gauge theories and Abelian gauge-theories with monopoles and leads to confinement 
when considered in Euclidean space-time . Therefore, we call the tensor structure 
multiplied by (1 — k) non-confining (nc) and the one multiplied by k confining (c). 
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With the correlators (|2.11| ) and ( |2.13| ) and the minimal surfaces shown in Fig. [I], 



the functions x P an d X 

d 2 k 



r NP 



Xnc Xc 



NP 



X 



47T 



a 8 (kt)iD'^(ki 



have the following form |1[ 

e ikj_(r lq -r 2q ) _|_ gifcx (r lq -r 2q ) 



_ e ik ±(ri q -r 2 q) _ e ik ± {f lq ~f 2q ) 



(2.18) 



NP 
Anc 



Tl 2 G 2 (l-K) f d 2 k 



24 



(2nf 



e ik ± (f lq -f 2q ) _|_ e ikt(r lq —r 2q ) 



_ e ikj_(riq-r2 9 ) _ e ik ± (r lq -r 2q ) 



(2.19) 



x N c p 



7T 2 G 2 K, 

24 



(ri ■ r 2 ) 



i ,i 



d 2 fc ± 

(2^ 



iD (2) (ki)e 



ik±(ri q +v 1 r 1 -r 2q -v 2 r 2 ) 



(2.20) 



where Z>p (2) (fci), D'i 2 \k]_) and D^l) are the Minkowskian correlation functions 
in transverse space obtained from (|2.11|) and ( 2.17 ), 



i 



±) 



iD?\kl) 



k 2 + m 2 G 

30 7T 2 

a(k 2 + a~ 2 

12 7T 2 

a (k 2 + a -2 



(2.21) 



f2.22) 



(2.23) 



The component x P describes the perturbative interaction of the quark and antiquark 
of one dipole with the quark and antiquark of the other dipole as evident from the 
r iq and dependence of (|2.18| ) and Fig. [L|b. The component Xnc nas the same 
structure as \ p and gives the non-perturbative interaction between the quarks and 
antiquarks of the dipoles. With the term quark used genuinely for quarks and 
antiquarks in the following, we refer to \ p and Xnc as quark- quark interactions. 

The component Xc given in ( |2.20| ) shows a different structure. Here, the inte- 
grations over v\ and v 2 sum non-perturbative interactions between the gluon field 
strengths connecting the quark and antiquark in each of the two dipoles. These 
connections are manifestations of the strings that confine the corresponding quark 
and antiquark in the dipole and are visualized in Fig. |l|b. Indeed, the SVMshows 
a flux tube between a static quark-antiquark pair in Euclidean space-time [|9| p7T| . 
Therefore, we understand the confining component x^ P as a string-string interac- 
tion. 
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The mixed contribution \nc Xc P that occurs in cross sections discussed in the 
next sections gives the non-perturbative interaction of the quark and antiquark of 
one dipole with the string of the other dipole, i.e., it represents the quark-string 
interaction. 

In our previous paper |TJ, we studied extensively the phenomenological perfor- 
mance of the model. With a strong and weak powerlike energy dependence ascribed 
respectively to \ p an d X NF '■> we have achieved a successful description of many dif- 
ferent observables for hadron-hadron, proton-photon and photon-photon reactions 
over a wide energy range. Therefore, we adopt from all given parameter val- 
ues. The phenomenological energy dependence, however, is not essential for the 
structural aspects discussed in this work. Observables obtained with the T-matrix 
element ([2.5|) agree with experimental data at a fixed cm. energy of a/sq ~ 20 GeV. 



In this paper, we focus on the analytic structure of the model and concentrate 
on the limit of small x-f unc ti° ns ) X P < 1 an d X NP < 1, in which the T-matrix 
element ( |2.5| ) simplifies to the sum of a perturbative (P) and non-perturbative (NP) 
component 



T{s ,t) = 2is d 2 b ± e 1 ^ d Zl d 2 n dz 2 d 2 r 2 |Vi(*i,ri)| 2 \ik{z 2 , r 2 )\ 



(x p ) 2 + (x Np y 



(2.24) 



We show in the next section that the perturbative component, (x ) , describes two- 
gluon exchange |19],|20| and that the non-perturbative component, (x wp ) 2 , represents 



the corresponding non-perturbative two-point interaction. The main motivation of 
this paper is to calculate the analytic structure of this non-perturbative interaction 
in momentum space. 



3 The Structure of Dipole-Dipole Scattering in 
Momentum Space 

In this section, we calculate the perturbative and non-perturbative QCD contribu- 
tion to the dipole-dipole {DD) scattering amplitude in momentum space. Beyond 
the known perturbative two-gluon-exchange interaction between the dipoles [[I9|,|2"0[| 
new insights into the non-perturbative structure of the scattering process are ob- 
tained within the SVM. The non-perturbative string-string interactions show a new 
structure different from the dipole factors of the perturbative contribution. 
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The total dipole-dipole cross section in the small- x limit is obtained from ( 2.24J) 
via the optical theorem 



1 

So 



lmT(s ,t = 0) 



2 d 2 b± dz i d 2 r 1 \dz 2 d 2 r 2 \^{z 1 ,r^)f\i) D ^{z 2 ,r 2 )\ A 



x- 



9 



NP i ^,NP\ 2 



(3.1) 



where the dipoles have fixed Z{ and \ fi\ values but are averaged over all orientations 

1 



(3.2) 



All integrations in ( |3.1| ) can be carried out analytically except the one over transverse 



momentum \k±\ that enters through the x-fmictions ( 2.18 )-( 2.20 ). The perturbative 
(P) and non-perturbative (NP) integrand of the resulting total dipole-dipole cross 
section 



d\k ± \ I F (s ,\k±\) + I NP (s ,\k±\) 



d\k ± \ r(s , \k ± \) + I^(s , \k ± \) + I^(s , \k ± \) 



(3.3) 
(3.4) 



show the following momentum-space structure 
I P (so,\kx\) 

2 (^Jl - e^|^) 2 (^ 2 |1 - e^ 2 |^ 2 ) 



2 -^\{A^a s (kl)f^bf\kl)] 2 



I?f(so,\k ± \) 



x 
2 1 



C>o,|*l|) 



9 2vr 

x 

2 1 



24 



{\- K)% bf\k])+^ i ib { - 2 \k 



-\ 2 



kl 



9 2vr 

x 



2 (Vkjl - e^lVbr) 2 (^ 2 |1 - e^ 2 |^ 2 ) 

n 2 

Ik 



24 



(3.5) 



(3.6) 



(3.7) 



2 (VbJ tan 2 0!(l - e lk ^)\^ Dl ) 2 (^ 2 | tan 2 2 (l - e^ 2 )|^ 2 ) 
where the brackets {ip Di \...\ip Di ) denote the averages 



(ip Di \Ai\ij; Di ) = / dz i d 2 r i \?p Di (zi,r i )\ 2 A i 



(3i 
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and the dipole orientation fa is defined as the angle between transverse momentum 
k± and dipole vector r*j. With (|3.2j ) the integration over the dipole orientations fa 
leads respectively to the Bessel function JoO&lII^dJ) an d the generalized hypergeo- 
metric function^ 1 F a (-l/2; 1/2, 1; — 







JbdfcLll^l) , (3-9) 



(Vytan^l-e^)!^) = ^- / ^ dfatan 2 fa(l - e ik ^) 





i i t 2 r 2 

= -l + iF 2 (-i;i 1;^) . (3.10) 

The important implications will be discussed in the next section. 

The integrand I p given in (|3.5| ) describes the known perturbative two-gluon| 
exchange [0,0 between the quarks and antiquarks of the two dipoles. The in- 



gredients of I are visualized for one combination of gluon exchanges in Fig. |2|a: 
the paired horizontal lines represent the dipole factors (1 — e tk± - ri ) that describe the 
phase difference between the quark and antiquark at separate transverse positions, 
the surrounding brackets indicate the average over the dipole orientations ( |3.8|) , the 
two curly lines illustrate the gluon propagator squared [iDp 2 ^] 2 , and the four vertices 
(dots) correspond to the strong coupling to the fourth power g 4 = (4na s ) 2 . 

The integrand I NP = I^ p + I^ p given in ( |3.6| ) and ( |3.7| ) describes the non- 
perturbative interactions: the quark-quark, string-string, and quark-string inter- 
actions identified by the appropriate correlation functions [iD[ ] 2 , [iD^] 2 , and 
[iD'^iD®]. These interact ions are illustrated in Figs. 0b, |2|c and [2]d, respectively. 
Analogous to the perturbative interaction in Fig. ||a, the dashed and solid zig- 
zag lines represent respectively the non-confining (1 — k) iD^ and the confining 
(n/kf) iD^ non-perturbative correlations, the shaded areas symbolize the strings, 
and the four vertices (squares) in each figure indicate the "non-perturbative coupling 
" to the fourth power g A NP := (7r 2 G 2 /24) 2 . 

The integrand I^ p describes the non-perturbative interactions between the quarks 
and antiquarks of the two dipoles and exhibits the same dipole factors (l—e tk±Tl ) that 



review of generalized hypergeometric functions can be found in |34(| . In the computer pro- 
gram Mathematica J3|| "HypergeometricPFQ [{ — 1/2} , {1/2,1}, denotes this function. 



The exact T-matrix element (2.5) important to describe scattering processes at ultra-high 



energies [jj] goes beyond two-gluon exchange due to the higher orders in the cosine expansion. 
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(a) (b) 




Figure 2: Perturbative and non-perturbative contributions to dipole-dipole scattering: (a) 
perturbative quark-quark interaction and non-perturbative (b) quark-quark, (c) string- 
string, and (d) quark-string interactions. The term quark is used genuinely for quarks 
and antiquarks. Only one diagram is shown for each type of interaction. Paired horizontal 
lines represent quark-antiquark dipoles, surrounding brackets indicate the averages over 
the dipole orientations ( |3.8[) , shaded areas visualize strings, curly lines denote exchanged 
perturbative gluons, and dashed and solid zig-zag lines symbolize respectively the non- 
perturbative non-confining and confining correlation functions. 



appear in the perturbative integrand ( |3.5|) . I^ q p contains three components: the non- 
confining component, the confining component, and their interference term. While 
the non-confining component visualized in Fig. |2|b has the same structure as the 
perturbative contribution, see also ( p.!8| ) and (|2.19|) , the confining component shown 
in its more general form in Fig. |^c comes from the interaction between the quarks 
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and antiquarks at the endpoints of the strings, which will be further discussed below 
Eq. ( [3. lip . The interference term describes the quark-string interaction as illustrated 
in Fig. |2|d. Note that it reduces entirely to an interaction between the quarks and 
antiquarks of the dipoles with the additional denominator l/kf_ generated by the 
integrations over the variables V\ and v 2 in the confining component \c (|2-20|) . 

The integrand iff describes the non-perturbative string-string interaction shown 
in Fig. |^c. The new angular dependencies in the string-string interaction, the mod- 
ified dipole factors tan 2 0j(l — e tk±rD ^) in (|3.7|) , are obtained as follows. The integra- 
tions over i>i, t> 2 , v[, and v' 2 in (xf P ) 2 produce the dipole factors (1 — e lk±n ) and the 
denominator l/((k±fi) 2 (k±r2) 2 ). This denominator multiplied with the additional 
factor if\f*2) 2 from (x^ P ) 2 , see ( P-20|) , gives the total angular dependence 

(rif 2 ) 2 r\r\ cos 2 (0i - 2 ) 



(4ri) 2 (4r 2 ) 2 (klr 2 cos 2 0!) (A; 2 r 2 cos 2 2 ) 

(cos 0i cos 02 + sin 0! sin 2 ) 2 



kj_ COS 2 0i COS 2 2 

= -j (l + 2 tan0x tan0 2 + tan 2 0! tan 2 2 ) . (3.11) 
± 

The first term in ( |3.11|) explains the interaction between the endpoints of the strings 
with the additional factor \fk\_ in the first contribution l^ p already mentioned 
above. The product of the second term in ( |3.11| ), 2tan0itan0 2 , with the dipole 
factors vanishes after the integration over the dipole orientations 0i and 2 . The 
third term in ( 3.11 ), tan 2 0i tan 2 2 , weights the different orientations of the strings 
and is characteristic for the string-string interaction ( |3.7|) in our model. Due to 
this factor, the string-string interaction differs significantly from the interaction 
between the quarks and antiquarks of the dipoles known from perturbative two- 
gluon exchange (fJ3|). 

In Fig. [| we show the perturbative integrand I p (solid line) and the non- 
perturbative integrand I NP (dashed line) of the total dipole-dipole cross section ( |3.4p 
as a function of transverse momentum \k±\ for various dipole sizes. The integrands 
have been calculated with the model parameters given in Sec. |2] that allow a good 
description of the experimental data of hadron and photon reactions at the cm. en- 
ergy of t/s~o ~ 20 GeV with the wave functions given in Appendix |A] and the exact 
T-matrix element ( |2.5| ). Evidently, the non-perturbative integrand I governs the 
low momenta and the perturbative integrand I p the high momenta exchanged in 
the dipole-dipole scattering. This behavior is a direct consequence of the correlation 
functions: The non-perturbative correlation functions (|2.22|) and ( |2.23| ) favor low 



momenta and suppress high momenta in comparison to the perturbative correlation 
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Figure 3: The perturbative integrand I p (solid line) and the non-perturbative integrand 
I NP (dashed line) of the total dipole-dipole cross section ( |3.4[ ) as a function of transverse 
momentum \k±\ for various dipole sizes \ f Dl \ and (r^l- The integrands are shown for the 
parameters given in Sec. [2] that allow a good description of the experimental data of hadron 
and photon reactions at ^/sq 20 GeV with the wave functions given in Appendix [A| 
and the exact T-matrix element (|2.5| ). The oscillations of the perturbative integrand Jr 
originate from the Bessel function Jo(| /£>»!)• 
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function fl2.ll]) . For fixed dipole sizes \r Dl \ and |fb 2 | the absolute values of the pertur- 
bative and non-perturbative integrand are controlled respectively by the parameters 
(m g , M) and (G%, a, k) given in Sec. ||. With decreasing dipole sizes the ratio 
jp jjNP j ncreases: f or large dipole sizes \ f Dl \ = \r D2 \ = lfm, the non-perturbative 
contribution gives the main contribution to the total dipole-dipole cross section; for 

small dipole S1ZGS I | — I ^i?2 

| = 0.1 fm, the perturbative contribution dominates. 

As can be seen from the analytic results (|3.5|)-( |3~T0l) , the total dipole-dipole 
cross section ( |3.4| ) or the forward scattering amplitude T(s , t = 0) does not depend 
on the longitudinal quark momentum fractions z Di of the dipoles. For t = the 
parameter z Di disappears upon the integration over Zi since only \if; Di (zi, f^) | 2 ( |3.2| ) 
depends on Zi. 

The structure presented for dipole-dipole scattering remains in reactions involv- 
ing hadrons and photons: the hadronic and photonic total cross sections are ob- 
tained from the total dipole-dipole cross section (|3~Tp by replacing \ip Di (zi, r*i)| 2 given 
in fl3.2p with the hadron and photon wave functions given in Appendix [A[ As the 
total dipole-dipole cross section, the total hadronic and photonic cross sections are 
independent of the parameters which control the Zj - distribution in the wave func- 
tions due to the normalization of the Z{ - distributions. The independence of the 
total hadronic cross section on the widths Azh can be seen immediately with the 
Gaussian hadron wave functions ( |A.1| ) 

Ml-e^l^h) = l-e~^ s ", (3.12) 
(^tan^(l-e^)l^) = -1 + + ^\k x \S h Erf ( 1^ ) (3.13) 



with the error function Erf(z) = J* dt exp(— t 2 ). These analytical results con- 

firm the Azh and z D - independence of the total dipole-proton cross section assumed 
in phenomenological models However, the non- forward hadronic scatter- 

ing amplitude T(so,t ^ 0) depends on the parameter Az^ as shown explicitly in 
Appendix [B[ Thus, the differential elastic cross section da el /dt(s,t) ( p.6| ) and its 
logarithmic slope B(s,t) (p.7|) are Azu - dependent. In fact, this Azh - dependence 
is essential for the agreement with experimental data (l[|. 

The | k± | - dependence of the perturbative and non-perturbative integrand for 
hadron-hadron, hadron-photon, and photon-photon cross sections at high photon 
virtualities is similar to the one of the perturbative and non-perturbative integrand of 
the dipole-dipole cross section shown in Fig. |3] for (\f Dl \ = \f D2 \ = 1 fm), (1^ | = 1 fm, 
|fb 2 | = 0.1 fm), and (\f Dl \ = \f D2 \ = 0.1 fm), respectively Of course, the absolute 
values differ and the oscillations of the perturbative integrand caused by the Bessel 
functions disappear. 
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4 Decomposition of the QCD String into Dipoles 
and the Unintegrated Gluon Distribution 



The unintegrated gluon distribution of hadronsf] !Fh( basic, universal quan- 

tity convenient for the computation of many scattering observables at small x. It 
is the central object in the BFKL [23| and CCFM [24"] evolution equations. Upon 
integration over the transverse gluon momentum \k±\ it leads to the conventional 
gluon distribution xGh(x,Q 2 ) used in the DGLAP evolution equation |2"5[] . The 
unintegrated gluon distribution is crucial to describe processes in which transverse 
momenta are explicitly exposed such as dijet |H] or vector meson |22J production 
at HERA. Its explicit \k±\ dependence is particularly suited to study the interplay 
between soft and hard physics. In this section, an exact representation of the string 
as a collection of stringless quark-antiquark dipoles is presented that allows to ex- 
tract the perturbative and non-perturbative contributions to J r h(x,k 2 ) from our 
total dipole-hadron cross section via \k±\ -factorization. 



We calculate the unintegrated gluon distribution J-"h{x,k^) with the T-matrix 
element in the limit of small ^-functions ( |2.24|) . Motivated by the successful de- 
scription of many data in our recent work [jXJ] we give a strong energy dependence 
to the perturbative contribution \ P an d a weak one to the non-perturbative contri- 
bution x NP -i 



(x 



P\2 



(X 



P\2 ( ^0 



X 



(x Np r - (\ v/, ) 2/ 



X 



(4-1) 



where the values of the exponents e p = 0.73 and e NP = 0.125 are adopted from |T]] 
and x = 2.4 • 10~ 3 is adjusted to reproduce at Q 2 = 1 GeV 2 the integrated gluon 
distribution of the proton xG p (x, Q 2 ) extracted from the HERA data [37 
small- x limit ( |2.24 ) considered here is applicable only for x = Q 2 /s > 10 
higher energies (lower Bjorken-x) the full T-matrix element (|2.5| ) has to be used in 
order to satisfy unitarity constraints of the S'-matrix M, RH] . 



The 
4 . At 



The total dipole-hadron cross section o~Dh(x, \r D \) is obtained from the total 
dipole-dipole cross section ( |3.1| ) by replacing \i/; D2 (z2, r 2 )| 2 with a squared hadron 
wave function \ip h {z2, Tq)\ 2 . Accordingly, the x-dependent total dipole-hadron cross 

9 The word hadron and the subscript h is used genuinely for hadrons and photons in this section. 
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section reads 



0Dh(x, \r D \) = / dkf 



(4.2) 



x 



+ 



+ 



(4vr« s (^)) 2 
n 2 G 2 



iDf\kl) ( 1 - M\k ± \\r D \) ) m\l - e^\^ h ) \ ( ^ 



X 



X 



24 



K 



iD^(kl) + (l-K)iD'^(kl 



K 

~kl 



iD^(kl) 



-1+ iF 2 (- , , 



1- Jo(NI»b|)J (^hll-e^l^) 
(Vhltau^l-e^)]^) 



1 1 l;^))(^|tan 2 2 (l-e^)|^)[(|) 



with the Bessel function Jo(|&l||^d|) an d the generalized hypergeometric function 
iF 2 (— 1/2; 1/2, 1; — fcjV^/4) derived in the previous section. 

For dipole sizes |r^| — > 0, the perturbative contribution to (JDh{ x ) I'd I) is known 
to vanish quadratically with decreasing dipole size, o~Dh(x, \r D \) oc rj. This behavior 
reflects the weak absorption of a small color-singlet dipole in the hadron and is 
known as color transparency. It can be seen immediately from ( |4.2| ) as 



(l- Jo(\k±\\r D \) 



for \r D 



and finite \k± 



(4.3) 



The non-perturbative contribution to crrjh(x, \f D \) gives color transparency as well 
since the generalized hypergeometric function behaves as 



1+ F( 11 1--^ 



k]_ r D 



for |fx»| ^0 and finite |/cl 



(4.4) 



For large dipole sizes, |fi,|>lfm, the perturbative contribution to cr Dh (x,\f D \) 
describing interactions of the quark and antiquark of the dipole with the hadron 
saturates since 

(l - J (|£l||'d|)J ~ 1 for large \k ± \\r D \ • (4.5) 

In contrast, the non-perturbative contribution to crDh{x, \v D \) increases linearly with 
increasing dipole size, aoh(x, \f D \) oc \r D \. This linear increase is generated by the 
interaction of the string of the dipole with the hadron: The string elongates linearly 
with the dipole size \f D \ and, thus, has a linearly increasing geometric cross section 
with the hadron. Indeed, this feature of the string can be seen analytically since 

f-l + lJ F 2 (-i;i oc |4||r fl | for large |4||r fl | . (4.6) 
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When considered in Euclidean space-time, the same string gives also the linear 



confining potential between a static quark and antiquark at large qq separations [TI| 
TTJ . Obviously, the behavior of the total dipole-hadron cross section is related to the 
confining potential. Furthermore, as we are working in the quenched approximation, 
there is no string breaking through dynamical quark- antiquark production at large 
dipole sizes. String breaking is expected to stop the linear increase of the total 
dipole-hadron cross section at dipole sizes of \r D \ ^ 1 fm analogous to the saturation of 
the static qq potential seen for large qq separations on the lattice in full QCD |39|, |40fl . 



In Fig. |], we show the perturbative (solid line) and non-perturbative (dashed 
line) contributions to the total dipole-proton cross section (Td p (x, \r D \) as a function 
of the dipole size \r D \ at x = xq = 2.4 • 10~ 3 , where the perturbative contribution 
is multiplied by a factor of 10. The dipole-proton cross section is computed with 
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Figure 4: Perturbative (solid line) and non-perturbative (dashed line) contributions to 
the total dipole-proton cross section ao P (x, \f D \) as a function of the dipole size \f D \ at 
x = xq = 2.4 • 10~ 3 . with the perturbative contribution multiplied by a factor of 10. The 
perturbative contribution shows color transparency at small \ f D \ and saturates at large 
\f D \. The non-perturbative contribution shows also color transparency at small \ f D \ but 
increases linearly with increasing \r D \. 
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the simple Gaussian proton wave function ( |A.l ) and illustrates the general features 
discussed above: The perturbative contribution shows color transparency at small 
\r D \ and saturates at large \ f D \. The non-perturbative contribution shows also color 
transparency at small \f D \ but increases linearly with increasing \f D \. 



Our result Q4.2| ) shows that the \k±\ - dependence of the hadron constituents fac- 
torizes from the rest of the process in both the perturbative and the non-perturbative 
contribution to the total dipole-hadron cross section. This factorization - known 
in perturbative QCD as \kx\ -factorization [^TJ - allows to define the unintegrated 



gluon distribution Th{x, k 2 ) as follows fi"2| pi 



47r 2 r 2 r [1 - Jo(\k±\\r D \) 

<?Dh(x, \r D \) = — jdkl r a s (k%)T h {x,kf 

6 J {\k±\\r D \) 2 



(4.7) 



For small dipole sizes \ f D \, this equation together with ( |4.3| ) and the integrated gluon 
distribution 

rQ 2 

xG h (x,Q 2 ) = / dklF h {x,kl) (4.8) 
Jo 



leads to the widely used perturbative QCD relation ES 



<7Dh(x, \r D \) = — ^ [a s (Q 2 )xG h {x,Q 2 )] 



■-c/rg 



2 ; 



(4.9) 



where c ~ 10 is estimated from the properties of the Bessel function Jq(|A;j_| \r D \) ^ . 



To extract the unintegrated gluon distribution !Fh{z, k±), we compare ( |4.7p with 
% using the following mathematical identity 



F( 11 1 . ~k± r D- 



dZ-[l-J (\k ± \\r D \Z) 



(4.10) 



As discussed in the previous section, the lhs of ( |4.10p results from the string averaged 
over all orientations Q3.10|) . Thus, the string confining the quark-antiquark dipole 
of length \r D \ can be represented as an integral over stringless dipoles of sizes £\r D \ 
with < C, < 1 and a dipole number density of n(£) = l/£ 2 . As visualized in 
Fig. the string-hadron scattering process reduces to an incoherent superposition 
of stringless dipole-hadron scattering processes with dipole sizes < £\r D \ < \r D \ 
and dipole number density n(£) = l/^ 2 . The decomposition of the string into many 
smaller stringless dipoles via ([4.1 0Q behaves similar to the wave function of a qq 
onium state in the large N c limit |47, 48|: The numerous gluons emitted inside the 



onium state can be considered as many qq dipoles in the large N c limit. 
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Figure 5: The string of length \f D \ is made up of stringless dipoles of size £\r D \ with 
< £ < 1 and dipole number density n(£) = l/£ 2 . The string-hadron scattering process 
reduces to an incoherent superposition of stringless dipole-hadron scattering processes. 



Inserting ( |4.10|) into ( |4.2|) and rescaling the momentum variable |A;_|| = £\k±\ 
the string-hadron (sh) contribution to the total dipole-hadron cross section ( f4~2] 
becomes 



°Dh(x, \r D \) 



x 



7T 2 G 2 

24 



9 4vr 

.2 pi 



dfcjMl- J (\kL\\r D \ 



(4.11) 



k' 



2 1 2 



(ip h \ tan 2 2 (l 



>l'v-)M7 



The dipole factor (1 — •/o(|&j_||^d|) indicates that the string-hadron interaction has 
been rewritten into a stringless dipole-hadron interaction. The string confining the 
dipole has been shifted into the hadron. Comparing (|4.7|) with our new expression 
for the total dipole-hadron cross section, 



^Dh(x,\f D \) = I dkl(l - J (\k ± \\f D \) 



x 



+ 



7T G 2 



24 

2 /•! 







:i-k! 



'^ (2) (f) 



iD^(kl) 



X 
l 2 



(i/j h \l-e lk ^\?p h ) 



(^|tan 2 2 (l-e^/^)|^)^- 



(4.12) 
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one obtains the unintegrated gluon distribution 

k'i 



+ 



+ 



K 



7T 2 G, 

24 

2 /■! 



2 fr 



(i-«) ^; (2) (^) 



+ 



(A\l-e lk ^ h ) 







i£>( 2 )( 



(^|tan 2 02(l-e^/^ 2 )|^) 



;r 



(4.13) 



This result shows explicitly the microscopic structure of the perturbative and non- 
perturbative contribution to T^ix, k±)- It is valid for any hadron wave function. To 
get numerical values for the unintegrated gluon distribution jF h (x, fcf), the hadron 
wave functions must be specified. 



Numerical Results for the Unintegrated Gluon 
Distribution in Hadrons and Photons 



In this section, we present the unintegrated gluon distribution ( |4.13| ) for protons, 



pions, kaons, and photons computed with the Gaussian hadron wave function ( A.l 
and the perturbatively derived photon wave functions ( JA.2| ) and ( |A.3[ ). To account 
for non-perturbative effects at low photon virtuality Q 2 in the photon wave functions, 
quark masses rrif(Q 2 ) are used that interpolate between the current quarks at large 
Q 2 and the constituent quarks at small Q 2 [|15| as discussed in Appendix [A]. 



The unintegrated gluon distribution of the proton J-p(x, kf_) is shown as a function 
of transverse momentum \k±\ at x — 10 _1 , 10~ 2 , 10 -3 , and 10~ 4 in Fig. |6|. Figure ^ 
illustrates the interplay of the perturbative (solid line) and non-perturbative (dashed 
line) contributions to \k±\J 7 p (x,kf) as a function of transverse momentum \k±\ for 
the same values of x. At small momenta \k±\, the unintegrated gluon distribution 
is dominated by the non-perturbative contribution and behaves as This be- 

havior reflects the linear increase of the total dipole-proton cross section at large 



dipole sizes. In contrast, the saturation model of Golec-Biernat and Wiisthoff | 44| 
shows the behavior J-^ BW (x, kf) cx k\ for small momenta. With increasing the 
non-perturbative contribution to J-j,(x, k±) decreases rapidly which results from the 
strong suppression of large momenta by the non-perturbative correlation functions 
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iD[ {2 \kl) and iD^(kl) given in ( gjg ) and ( |2T23|) . For |£j > 1 GeV, the perturba- 
tive contribution dominates the unintegrated gluon distribution. It drops as l/k]_ 
in accordance with the perturbative correlation function iDp 2 \k]_) given in fl2.21|) . 
This perturbative QCD result is not reproduced by the phenomenological model of 
Golec-Biernat and Wiisthoff [44] which predicts a Gaussian decrease of T p (x, kj) 
with increasing \k±\. 
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Figure 6: The unintegrated gluon distribution of the proton T p {x,kJ) as a function of 
transverse momentum \k±\ at Bjorken-x values of 10 _1 , 10 -2 , 10 -3 and 10 -4 . 



The x- dependence of J-p(x,k1) can be seen in Figs. || and |7|: With decreasing 
x, the perturbative contribution increases much stronger than the non-perturbative 
contribution which results from the energy exponents e p 3> in (|4.1|) necessary 
to describe the data within our model Moreover, the perturbative contribution 
extends into the small - \ k±\ region as x decreases. Indeed, the soft-hard transition 
point moves towards smaller momenta with decreasing x as shown in Fig. [7[ Such 
a hard-to-soft diffusion is observed also in |E] where the unintegrated gluon dis- 
tribution has been parametrized to reproduce the experimental data for the proton 
structure function F2(x,Q 2 ) at small x. The opposite behavior is obtained in the 
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Figure 7: The unintegrated gluon distribution of the proton T p {x,kJ) times the tranverse 
momentum as a function of \k±\ at Bjorken-x values of 10 , 10~ 2 , 10~ 3 and 10~ 4 . 



approach of the color glass condensate fl50fl : With decreasing x, gluons are pro- 
duced predominantely in the high -\k±\ region of lower density and weaker repulsive 
interactions. 

In Fig. H, the unintegrated gluon distribution of the proton, pion, and kaon 
!Fh(x,kf} times the transverse momentum \k±\ is shown as a function of \k±\ at 
x = 10 -3 . The hadrons are characterized by different values for Azh and Sh in the 
hadron wave function (|A.1| ). However, J~h(x, kf) depends only on Sh- Due to the 
normalization of the hadron wave functions, A4 disappears upon the integration 
over Zi as can be seen directly from (|3.12| ) and ( |3.13|) . At small momenta, J~h(x, k]_) oc 
Si is found from (|3.12|) , (|3.13[) , and ( |4.13|) . It becomes visible in Fig. || for the 
chosen hadron extensions: S p = 0.86 fm, S n = 0.607 fm, and Sk = 0.55 fm. At 
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Figure 8: The unintegrated gluon distribution of the proton, pion, and kaon J r j l {x,k'^ i 
times the transverse momentum as a function of at Bj or ken- variable x = 10~ 3 . 



large momenta where the perturbative contribution dominates, the dependence on 
Sh vanishes as can be seen from ( |3.12| ) and the unintegrated gluon distributions of 
protons, pions, and kaons become identical. Of course, this behavior results from 
the wave function normalization being identical for protons, pions, and kaons with 
two valence constituents which are the quark and antiquark in the pion and kaon 
and the quark and diquark in the proton. At large \k±\, i.e., high resolution, the 
realistic description of protons as three-quark systems becomes necessary. In fact, 
the three-quark description of protons leads to a different perturbative contribution 
in (|4.13| ): the quark counting factors of 2 (appropriate for mesons) in the square 
brackets in Q3.5D are substituted by factors of 3 (appropriate for baryons) [19|,|20|]. At 
other values of x, the unintegrated gluon distributions of protons, pions and kaons 
show the same features. 

Photons are particularly interesting because the transverse size of the quark- 
antiquark dipole into which a photon fluctuates is controlled by the photon virtual- 
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Figure 9: The unintegrated gluon distribution of the transverse polarized photon 
T^* (x, k 2 , Q 2 ) times the tranverse momentum \k±\ as a function of \k±\ at photon vir- 
tualities of Q 2 = 0, 1, 10, and 100 GeV 2 and Bj or ken- variable x = 10~ 3 . 

ity Q 2 (cf. Appendix [A]) 



1 Q 2 + 4ml(Q 2 ) ' y - J 

where m u (Q 2 ) is the running w-quark mass given in Appendix [A]. In Figs. [5] and [II], 
the unintegrated gluon distribution of transverse (T) and longitudinal (L) photons 
J-^ (x, k]_, Q 2 ) times the transverse momentum is shown as a function of \k±\ 
for various photon virtualities Q 2 at x = 10 -3 . With increasing Q 2 , i.e., decreasing 
"photon size" |r 7 * |, the ratio of the perturbative to the non-perturbative contri- 
bution to J^,* (x, k 2 , Q 2 ) increases. This behavior has already been discussed in 
Sec. |3] on the level of dipole-dipole scattering, where decreasing dipole sizes increase 
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Figure 10: The unintegrated gluon distribution of the longitudinally polarized photon 
^^(x, k 2 , Q 2 ) times the tranverse momentum \k±\ as a function of | A£j_ | at photon virtual- 
ities of Q 2 = 0.01, 1, 10, and 100 GeV 2 and Bj or ken- variable x = IO" 3 . 



the ratio of the perturbative to non-perturbative contribution, see Fig. |^. Due to 
the different Q 2 - dependence in the transverse and longitudinally polarized pho- 
ton wave functions, (x, k 2 , Q 2 ) decreases continously with increasing Q 2 while 
T-^ix, k 2 , Q 2 ) increases for Q 2 ^l GeV 2 and decreases for Q 2 > 1 GeV 2 . 

In Fig. the unintegrated gluon distributions of hadrons and photons dis- 
cussed above are compared. The unintegrated gluon distribution of real (Q 2 = 0) 
photons (\fj \ ~ Sh) is suppressed by a factor of order a in comparison to the one 
of the hadrons otherwise its shape is very similar. The suppression factor comes 
from the photon- dipole transition described by the photon wave functions given in 
Appendix [X]. The ratio of the perturbative to the non-perturbative contribution 
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Figure 11: The unintegrated gluon distribution of the proton, pion, kaon, and transverse 
and longitudinally polarized photon Th,^* L {x, kf) times the transverse momentum \k±\ as 

a function of \k±\ at Bj or ken- variable x = 10 -3 . Results for transverse polarized photons 
are shown for photon virtualities of Q 2 = and 10 GeV 2 and results for longitudinally 
polarized photons for photon virtualities of Q 2 = 0.01 and 10 GeV 2 . 



of unintegrated gluon distributions increases as one goes from hadrons to photons 
with high virtuality Q 2 . Since the wave functions of protons, pions, and kaons are 
normalized to the same value, the unintegrated gluon distributions of these hadrons 
are, as mentioned above, identical at large \k±\ and do not depend on the hadron 
size. In contrast, the Q 2 - dependence of the photon wave functions leads to Q 2 - 
dependent, i.e., "photon size" - dependent, unintegrated gluon distributions at large 
\k±\. With increasing \k±\, the unintegrated gluon distributions of hadrons and pho- 
tons become parallel in line with the vanishing dependence on the specific form of 
the wave function. 



In Fig. [12], the integrated gluon distribution of the proton xG p (x, Q 2 ) is shown 
as a function of Bjorken-a; at photon virtualities of Q 2 = 1, 5, and 20 GeV 2 . Recall 
that the parameter xq = 2.4 • 10 -3 has been adjusted in the previous section such 
that the experimental data of xG p (x,Q 2 ) at Q 2 = 1 GeV [37 are reproduced. For 
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Figure 12: The gluon distribution of the proton xG p (x,Q 2 ) as a function of Bjorken-x 
at photon virtualities of Q 2 = 1, 5, and 20GeV 2 . 

x > 10~ 3 , xG p (x, Q 2 ) is mainly determined by non-perturbative physics as can be 
seen from Figs. [| and [7]. Perturbative physics becomes relevant for x ^ 10~ 3 and 
generates the steep increase of xG p (x, Q 2 ) with decreasing x at fixed Q 2 . Also 
the rise of xG p (x, Q 2 ) with increasing Q 2 at fixed x results from the perturbative 
contribution. For x <C 10~ 4 , we have shown explicitly in our recent work that 
multiple gluonic exchanges contained in the full T-matrix element ( p.5|) slow down 
the powerlike increase of xG p (x, Q 2 ) with decreasing x in accordance with S'-matrix 
unitarity constraints. 
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6 Comparison with Other Work 



In this section, we compare the unintegrated gluon distribution of the proton ex- 
tracted from our loop-loop correlation model (LLCM) with those obtained from the 
saturation model of Golec-Biernat and Wiisthoff (GBW) |44|], the derivative of the 
Gliick, Reya, and Vogt (GRV) parametrization of xG p (x, Q 2 ) [j5lH , and the approach 
of Ivanov and Nikolaev (IN) 



In the approach of Golec-Biernat and Wiisthoff [P], the unintegrated gluon 
distribution of the proton is extracted from the total dipole-proton cross section by 
inverting ( f4.7|) , 



J- p(x, kj_ 



3cr 
47r 2 a. 



Rl{x) 



k]_ exp (— Rl(x) kl) with 



(6.1; 



where the parameter a = 29.12 mb, a s = 0.2, Q = 1 GeV, A = 0.277, and x = 
0.41- 10~ 4 are obtained from a fit to the proton structure function F 2 (x, Q 2 ) including 
charm quarks in the photon wave function. Note, however that the GBW approach 
uses the dipole-proton cross section of the GBW model on the lhs of (|4.7| ) which 
implies multiple gluon exchanges while the rhs of ( f4.7| ) describes only two-gluon 
exchange as discussed in Sec. §. Moreover, as demonstrated below, the large k 2 - 
behavior of the unintegrated gluon distribution fl6.1| ) deviates significantly from the 
DGLAP results. This mismatch motivated the recent modifications of the saturation 



model [52 1 



The unintegrated gluon distribution of the proton is also computed from the 
integrated gluon distribution xG p (x, Q 2 ) by inverting (|4.8|) 



dxG p (x, Q 2 



dQ 2 



(6.2) 



Here, we use for the integrated gluon distribution xG p (x, Q 2 ) the leading order (LO) 



parametrization of Gliick, Reya, and Vogt (GRV) ||51|| , which covers the kinematic 
region 10" 9 < x < 1 and 0.8 GeV 2 < Q 2 < 10 6 GeV 2 . 



Ivanov and Nikolaev j|9| have constructed a two-component (soft + hard) ansatz 
for the unintegrated gluon distribution of the proton 



k 2 + tf 



k 2 + k 2 



(6.3) 
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with the soft and hard component 



•F 3oft(^J. 



Osoft Cp N c 



a s (kj 



71 



1 + 



3_kl 
A 2 



(6.4) 



* 1 J ^(05) 



9(Q'-*x)+J r p«(i,^)9(*x-Q') (6.6) 



pfH^i) nas the same form as flfcj.4| ) with the parameters a pt and /i pt in- 
stead of a so f t and yU so ft and J r pt (x,k^) is the derivative of the integrated gluon dis- 



where Tl\ (k 



tribution ( fT2|) . In the IN approach, the running coupling is given by a s {kj_) = 
min{0.82, (4tt)/(/? log^/A^])}. With the GRV-parametrization |T| for the 
integrated gluon distribution, the structure function of the proton F 2 (x,Q 2 ) has 
been described successfully using the following parameters: k 2 = 3 GeV 2 , k\ = 
(1 + 0.0018 log(l/x) 4 ) 5 , a soft = 2, a pt = 1, // soft = 0.1 GeV, /i pt = 0.75 GeV, 
Ql = 0.895 GeV 2 , (3 = 9, and A QCD = 0.2 GeV. 



In Fig. [13|, we show the LLCM, GBW, GRV, and IN results for the unintegrated 
gluon distribution of the proton T p {x, kf_) as a function of transverse momentum 
squared k\ for x = 10 _1 , 10~ 2 , 10~ 3 , and 10~ 4 . At small transverse momenta, 
k]_^j 0.1 GeV 2 , our model gives the largest values for T v {x^ kj_). As mentioned in 
the previous section, our LLCM unintegrated gluon distribution increases as 1/ 
with decreasing k]_ as a consequence of the linear increase of the total dipole-proton 
cross section at large dipole sizes. In contrast, for k\ — > 0, the unintegrated gluon 
distribution of GBW decreases as kj_ and the one of IN as k^_. In the perturbative re- 
gion, k]_^l GeV 2 , the unintegrated gluon distribution of the LLCM becomes smaller 
than the one of GRV and IN but is still larger than the one of GBW. Moreover, the 
LLCM, GRV, and IN unintegrated gluon distributions become parallel for x ^ 10~ 2 
and drop as l/k]_ for large k\. This perturbative QCD behavior is not reproduced 
by the GBW unintegrated gluon distribution which decreases exponentially with 
increasing k]_. 



The x-dependence of the LLCM, GBW, GRV, and IN unintegrated gluon dis- 
tributions T v [x^k^) is shown for transverse momenta squared k\ = 0.1, 0.5, 10, 
and 50 GeV 2 in Fig. [14]. In the non-perturbative region, k\ = 0.1 GeV 2 and 
k]_ = 0.5 GeV 2 , the LLCM, GBW, and IN unintegrated gluon distributions show a 
weak X- dependence. In the perturbative region, k\ = 10 GeV 2 and k\ = 50 GeV 2 , 
the x-dependence of the unintegrated gluon distributions becomes stronger. The 
LLCM, GRV, and IN unintegrated gluon distributions show nearly the same rise 
with decreasing x. In contrast, the GBW unintegrated gluon distribution increases 
much faster as x decreases. 
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Figure 13: The unintegrated gluon distribution of the proton J 7 p ( function of 

transverse momentum squared kj_ at Bjorken-x values of 10 _1 , 10~ 2 , 10 -3 , and 10 -4 . 
The different curves are obtained from our loop- loop correlation model (LLCM), the 
Golec-Biernat and Wiisthoff (GBW) model [44|, the derivative of the Gliick, Reya, and 
Vogt (GRV) parametrization of xG p (x,Q 2 ) 51], and the Ivanov and Nikolaev (IN) ap- 
proach 
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Figure 14: The unintegrated gluon distribution of the proton J T p ( function of 

Bjorken-x at transverse momenta squared kj_ = 0.1, 0.5, 10, and 50 GeV 2 . The different 
curves are obtained from our loop- loop correlation model (LLCM), the Golec-Biernat 
and Wiisthoff (GBW) model [13], the derivative of the Gluck, Reya, and Vogt (GRV) 
parametrization of xG p (x, Q 2 ) 5^], and the Ivanov and Nikolaev (IN) approach f49| . Note 
that the GRV parametrization is only available for k]_ > 0.8 GeV 2 . Moreover, the GBW 
result for F p (x,kj_) is below 10 -5 for k\ = 50 GeV 2 . 
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In addition to the one-scale unintegrated gluon distributions, J-{x, kj_), discussed 



in this work, there exist also two-scale unintegrated gluon distributions, T{x, kj_, /jr). 
In the CCFM evolution equation ||24|| , this additional scale /i 2 is related to the max- 
imum angle allowed in the gluon emission. Two-scale unintegrated gluon distribu- 
tions are obtained in the approach of Bliimlein [|53j| , Jung and Salam [54], Kimber, 
Martin, and Ryskin [55], and in the linked dipole chain (LDC) model [0,0. A 
comparison of their results can be found in |57|,|5!| where also the one-scale unin- 
tegrated gluon distributions of Kwiecinski, Martin and Stasto [59] , and Ryskin and 
Shabelski ISO) are discussed. 



7 Conclusions 



We have investigated perturbative and non-perturbative QCD contributions to high- 
energy scattering in momentum space within the recently presented loop-loop cor- 
relation model [0]. The perturbative contribution has been described by two-gluon 
exchange and the non-perturbative contribution by the stochastic vacuum model 
of QCD which leads to confinement of the quark and antiquark in a dipole via a 
string of color fields. This QCD string gives important non-perturbative contribu- 
tions to high-energy reactions and manifests itself in the linear increase of the total 
dipole-hadron cross section for large dipole sizes. A new structure different from 
perturbative two-gluon exchange has been found in the non-perturbative string- 
string interaction. A decomposition of the confining string into dipoles has allowed 
us to extract the microscopic structure of the unintegrated gluon distribution of 
hadrons and photons from the dipole-hadron and dipole-photon cross section via 
\k±\- factorization. 

The minimal surfaces used in our loop-loop correlation model [p]] have allowed us 
to show for the first time the QCD structure of the non-perturbative contribution 
to the dipole-dipole scattering amplitude in momentum space. This contribution 
has two parts: The first part describes the non-perturbative interaction between 
the quarks and antiquarks of the two dipoles and has the structure known from 
perturbative two-gluon exchange |n| |20| . The second part describes the interaction 
between the strings of the two dipoles and has a new structure originating from 
the geometry of the strings. The string contribution, however, is negligible in the 
scattering of two small dipoles which is governed by perturbative physics. 

The contribution of the confining string to the total dipole-hadron cross section 
&Dh(x, \r D \) has been studied. For small dipole sizes, |fb| — > 0, the string contribution 
shows color transparency, <JDh{x, \r D \) known for the perturbative contribu- 

tion. For large dipole sizes, |r^| ^ 0.5 fm, the non-perturbative contribution increases 
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linearly with increasing dipole size, <JDh{x, \r D \) oc \f D \, in contrast to the perturba- 
tive contribution which gives an \ f D \ - independent dipole-hadron cross section. This 
linear increase is generated by the interaction of the string of the dipole with the 
hadron: the longer the string, the larger the geometric cross section with the hadron. 
String breaking, expected to stop the linear increase for \ f D \ <; 1 fm, is not obtained in 
our model working in the quenched approximation. In our model, \k± \ - factorization 
has been found to be valid also for the non-perturbative contribution to (TDh(x, |^d|)- 

A very interesting feature of the string confining the quark and antiquark in the 
dipole has been found: A string of length \r D \ can be represented as an integral over 
stringless dipoles of sizes £\r D \ with < £ < 1 and dipole number density n(£) = 
l/£ 2 . A similar behavior has been observed for the perturbative wave function of a qq 
onium state in the large - N c limit where the numerous emitted gluons are considered 



as dipoles |47||48[. The decomposition of the string into stringless dipoles has allowed 
us to rewrite the string-hadron scattering process as an incoherent superposition of 
dipole-hadron scattering processes and to extract the unintegrated gluon distribution 
of hadrons and photons from our dipole-hadron and dipole-photon cross section via 
\k±\- factorization. 

We have shown explicitly the microscopic structure of the unintegrated gluon 
distribution in hadrons and photons. For small momenta \k±\, the unintegrated 
gluon distributions of protons, pions, and kaons T^x^kf) are dominated by non- 
perturbative physics and behave as S%/\k±\ where Sh denotes the hadron extension. 
The l/|fc±| - behavior reflects the linear increase of the total dipole-hadron cross sec- 
tion at large dipole sizes. The unintegrated gluon distributions of protons, pions, and 
kaons differ at small momenta because of their different extensions: S p = 0.86 fm, 
S n = 0.607 fm, and Sk = 0.55 fm. For large momenta, <^ 1 GeV, the uninte- 
grated gluon distributions of the hadrons are dominated by perturbative physics 
and show the - behavior induced by the gluon propagator. In the perturbative 
region of large momenta, the valence constituents are resolved and the dependence 
of the unintegrated gluon distribution on the hadron extension Sh vanishes. In con- 
trast, the unintegrated gluon distribution of photons depends on the "photon size" 
controlled by the photon virtuality Q 2 for large \k±\. For very large \k±\, the uninte- 
grated gluon distributions of hadrons and photons become parallel in line with the 
vanishing dependence on the specific form of the wave function. 

The x- dependence of the unintegrated gluon distribution of hadrons and photons 
has been introduced phenomenologically into our model. Motivated by the successful 
description of many experimental data in our recent work |T|, we have given a 
strong energy dependence to the perturbative contribution and a weak one to the 
non-perturbative contribution. Consequently, with decreasing x the perturbative 
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contribution increases much stronger than the non-perturbative contribution and 
extends into the small - \k±\ region. A similar hard-to-soft diffusion is observed also 
in the approach of Ivanov and Nikolaev ||9| while a soft-to-hard diffusion is obtained 



in the approach of the color glass condensate ||50|| . Considering the integrated gluon 
distribution of the proton xG p (x, Q 2 ), our non-perturbative contribution dominates 
for x ?t 10 -3 while the perturbative contribution becomes relevant for x ^ 10~ 3 and 
generates the steep increase of xG p (x, Q 2 ) with decreasing x at fixed Q 2 . Also 
the rise of xG p (x, Q 2 ) with increasing Q 2 at fixed x results from the strong energy 
dependence of the perturbative contribution. 

We have compared the unintegrated gluon distribution of the proton T p {x, k 2 ) 
extracted from our loop-loop correlation model (LLCM) with those obtained from 
the saturation model of Golec-Biernat and Wiistlioff (GBW) [[]] . the derivative 
of the Gluck, Reya, and Vogt (GRV) parametrization of xG p (x,Q 2 ) |5l|, and the 
approach of Ivanov and Nikolaev (IN) For k 2 — > 0, the unintegrated gluon 



distribution of GBW decreases as k 2 and the one of IN as fef in contrast to the 
1/ a/aJ- decrease found in our model. In the perturbative region, the LLCM, GRV, 
and IN unintegrated gluon distributions become parallel for x ^ 10 -2 and drop as 
\jk\ with increasing k\. This perturbative QCD behavior is not reproduced by the 
GBW unintegrated gluon distribution which decreases exponentially with increasing 
k\. The x-dependence of the considered unintegrated gluon distributions is weak in 
the non-perturbative region and becomes stronger as kj_ increases. 

The unintegrated gluon distributions of hadrons and photons, J r h(x,k 2 ), ex- 
tracted from our model can be used to compute, for example, dijet and vector 
meson production. So far, the unintegrated gluon distribution is obtained from the 
lowest order contribution to the T-matrix element and can only be used for x ?t 10~ 4 . 
For smaller values of x, multiple gluonic exchanges become important and necessary 
to respect S-matrix unitarity constraints jl]. Thus, a formalism must be developed 
that allows to extract an unintegrated gluon distribution which takes into account 
multiple gluonic exchanges. 
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A Hadron and Photon Wave Functions 



The light-cone wave functions ipi(zi,ri) provide the distribution of transverse size 
and orientation r*j and longitudinal quark momentum fraction z% to the light-like 
Wegner- Wilson loops VK[Cj] that represent the scattering color-dipoles. 



The Hadron Wave Function 

In this work, mesons and baryons are assumed to have a quark- ant iquark and 
quark-diquark valence structure, respectively. This allows us to model hadrons as 
color-dipoles. We use for the hadron wave function the phenomenological Gaussian 
Wirbel-Stech-Bauer ansatz 



where the constant Nh is fixed by the normalization of ( A.lj) to unity. The dif- 



ferent hadrons considered - - protons, pions, and kaons - - are characterized by 
different values for Azh and S^. The extension parameter is a fit parameter 
that should resemble approximately the electromagnetic radius of the correspond- 
ing hadron while Azh = w/(y2mf l ) PTJ * s fi xe d by the hadron mass and the 



value w = 0.35 — 0.5 GeV extracted from experimental data. We adopt the values 
of Az h and Sh from our previous work |JXJ] as these values allow a good description 
of many experimental data: Az p = 0.3 and S p = 0.86 fm for protons, Az^ = 2 and 
S n = 0.607fm for pions, and Azk = 0.57 and Sk = 0.55 fm for kaons. 



The Photon Wave Function 



The photon wave function ip-yfa, fi, Q ) describes the fluctuation of a photon with 
virtuality Q 2 into a quark-antiquark pair with longitudinal quark momentum frac- 
tion Z{ and spatial transverse size and orientation r*j. The computation of the cor- 
responding transition amplitude (qq(Zi,fi)\j(Q 2 )) can be performed conveniently in 
light-cone perturbation theory |HJ] and leads to the following squared wave functions 



for transverse (T) and longitudinally (L) polarized photons ||2| 
3 o 
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hM^,rl,Q 2 )| 2 = I^E^iW + a-^l^^/I^O+^/^/lrll)} (A.2) 

/ 



l^fe^Q 2 )! 2 = |^5> 2 {AQ 2 z 2 (l-z l ) 2 K 2 (e f m}, (A.3) 



2tt 2 

/ 
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where a is the fine-structure constant, ej is the electric charge of the quark with 
flavor /, and K and K\ are the modified Bessel functions (McDonald functions). 
In the above expressions, 

e) = Zi (l- Zi )Q 2 + m) (A.4) 

controlls the transverse size (-distribution) of the emerging dipole, \fi\ oc 1/e/, that 
depends on the quark flavor through the current quark mass mf. 



For small Q 2 , the perturbatively derived wave functions, ( |A.2| ) and (|A.3| ), are 



not appropriate since the resulting large color-dipoles, i.e., \fi\ oc 1/rrif ^> 1 fm, 
should encounter non-perturbative effects such as confinement or chiral symmetry 
breaking. To take these effects into account, we introduce Q 2 -dependent quark 
masses, rrif = m/(Q 2 ), that interpolate between the current quarks at large Q 2 and 
the constituent quarks at small Q 2 [jnj. In our recent work |], we adjusted the 
running quark masses as follows 

m M (Q 2 ) = 0.178 GeV (1 - ®{Q 2 u4 - Q 2 ) , (A.5) 

o 2 

m s (Q 2 ) = 0.121 GeV + 0.129 GeV (1 - 6(Q 2 - Q 2 ) , (A.6) 

Qt 

with the parameters Q 2 l d = 1.05 GeV 2 and Q 2 = 1.6 GeV 2 . For the charm quark 
mass, we used a fixed value of 

m c = 1.25 GeV . (A.7) 



B The Non-Forward T-Matrix Element 

In this appendix, we calculate the perturbative and non-confining contribution to the 
non-forward (t ^ 0) T-matrix element. We show explicitly that the non-forward T- 
matrix element depends on the parameters which control the Zi - distribution of the 
wave functions. These parameter are essential for a good description of differential 
elastic cross sections and the slope parameter as shown in our recent work 

The confining contribution to the non-forward (t ^ 0) T-matrix element is not 
presented. It is much more complicated since some of the integrations cannot be 
performed analytically. Nevertheless, we find numerically that it shows the same 
features concerning the Z{ - distributions as the perturbative and non-confining con- 
tributions. 
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The perturbative contribution to the non-forward T-matrix element in the small- 
X limit 



T p (s ,t) = ^° ld 2 b ± e^ ld Zl d 2 ri ldz 2 d 2 r 2 \M^un)\ 2 |^ 2 ,r 2 )| 2 { X P f , 
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reduces upon integration over the impact parameter to 



(B.1J 



T P (s ,t) 



32is ( 
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d 2 k ± a s (k 2 ) iDf\kl) aMk ± + ql) 2 ) iD'f> Uk ± + a± 



V(2) 



1 r l 



H x (zfql) - H x + k ± ) 2 ) H 2 (z 2 2 ql) - H 2 Uz 2 q\_ + k 



x / dz\ \ dz 2 
'0 Jo 



with 

#i(W± + 4)) := Jd 2 ri \Uzi,ri)\ 2 e^^ + ^ 
and \ipi(zi,fi)\ 2 denoting hadron or photon wave functions. 



(B.2) 
(B.3) 



The non-confining contribution to the non-forward T-matrix element in the 
small- x limit ( |2.24 ), 



2tSr 



T/rU,./-) = =^ ld 2 b ± e^ \dz x d 2 r x \dz 2 d 2 r 2 \M^,n)\ 2 \M^r 2 )\ 2 ( X ^) 2 



becomes analogously 



is ( n 2 G 2 (1 — k) 



d 2 k\ 



(B.4) 
(B.5) 



x / dz\ I dz 2 
'0 Jo 



Ht (zlql) - H x {(z iq l + k ± ) 2 j j \H 2 (zjql) - H 2 {(z 2 ql + k±) 
with Hi 2 defined in (|B.3|) . 



For t = —q]_ 7^ 0, both contributions ( |B.2| ) and ( |B.5[ ) depend on the shape of 
the Zj-distribution of the wave function, i.e., for the Gaussian hadron wave func- 
tion flA.lp , the contributions depend on the width Az h . This Az ft -dependence is 
transferred to the differential elastic cross section 



da 



el 



dt 



(s,t) 



16tts' 



■\T(sM 



and its local slope 



fl(.,«) = s lln 



da 



el 



dt 



(s,t) 



(B.6) 



(B.7) 
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At t = 0, the dependence on the shape of the ^-distribution of the wave functions 
disappears because of the normalization of the ^-distribution as can be seen imme- 
diately from flB.3|) . Therefore, in our model the total cross sections - related via the 
optical theorem to the forward (t = 0) T-matrix element - do not depend on the 
parameter that characterize the Zj-distribution of the wave function. 
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